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Abstract. The Lattice Hadron Physics Collaboration (LHPC) baryon spectroscopy effort is re- 
viewed. To date the LHPC has performed exploratory Lattice QCD calculations of the low-lying 
spectrum of Nucleon and Delta baryons. These calculations demonstrate the effectiveness of our 
method by obtaining the masses of an unprecedented number of excited states with definite quan- 
tum numbers. Future work of the project is outlined. 



A main goal of the LHPC is to determine the spectrum of excited baryons. This 
is achieved by extracting hadron masses from the exponential fall-off of two-point 
temporal correlation functions (0| @i{t) ^j(O) |0) , in which the operators G{ are composed 
of quark and gluon fields. These correlators can be expressed in terms of path integrals 
on a discrete space-time lattice [1] which are calculated numerically using Monte Carlo 
techniques [2]. 

The method developed by the LHPC to extract baryonic spectra is detailed in Refs. [3, 
4, 5]. This involves the construction of interpolating operators to create states with 
definite quantum numbers, the evaluation of the temporal correlators (0|^/(f)^/(0)|0), 
and the pruning of the final operator sets. 

If we have a set of N operators we may form the N xN (time-dependent) matrix 

Cij{t) = (O\0i(t)0j(O)\O). (1) 

The diagonalization of C~ 1//2 (fo)C(/)C~ 1//2 (fo), where to is a small non-zero reference 
time, produces correlators that may be fit to decaying exponentials to obtain Ej, the ener- 
gies of the lowest-lying states that can be interpolated by {^/}. If only zero momentum 
operators are used, the extracted energies of single particle states correspond to their 
masses. 
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FIGURE 1. The various types of extended baryon operators. Solid circles represent smeared quark 
fields, lines represent smeared link fields, and hollow circles the location of the reference site. 



Our baryon operators are composed of covariantly displaced quark fields as build- 
ing blocks combined in a gauge invariant way to interpolate both ground and excited 
baryonic states. Both the quark and link fields are smeared (replaced with a spatially 
localized weighted average) [4] to increase overlap with low-lying modes. Furthermore, 
they are designed to be spatially extended in order to interpolate states which are radially 
and orbitally excited. This spatial extension is achieved by gauge-covariantly displacing 
some (or all) of the three quarks in different directions from some reference site. The 
various types of extended baryon operators are shown in Fig. 1. 

Group theoretical projections are used to construct operators which transform irre- 
ducibly under the symmetries of the lattice. Recall that a baryonic state is denoted by its 
mass (m), spin and parity (J p ), and isospin (I). To construct states with definite isospin, 
we combine quark fields of different flavors to create an operator that transforms ir- 
reducibly under SU(2) isospin. Since the QCD Lagrangian is diagonal in flavor (and 
thus isospin) space, members of an iso-multiplet will be mass degenerate. Thus we may 
arbitrarily choose the operator with maximal third component of isospin 1$ as the inter- 
polating operator for a given iso-multiplet. 

In continuous space, J p labels an irreducible representation of the (double-valued) 
rotation group SU(2), which is a symmetry group of rotationally invariant Hamiltonians. 
However, on a cubic lattice the symmetry group is no longer that of continuous rotations 
but the subgroup of discrete lattice rotations. This group is denoted and is termed 
the double cubic point group. Because this group has a finite number of elements, 
it also has a finite number of irreducible representations, in contrast to the infinite 
number of possible J p values in continuous space. These representations are labeled 
Gig,Hg,G2g,Gi u ,H u , and G2 U for historical reasons. There also exists a procedure 
known as subduction to identify the states appearing in these irreducible representations 
with continuum J p values. For a complete review of the operator construction procedure, 
seeRef. [3]. 

After employing these considerations, we have operators which interpolate states 
with definite isospin, lattice spin, and parity. However, for a particular isospin (e.g. 
the Nucleons or Delta baryons) and lattice spin-parity value (such as Gi g or H g ), this 
operator set is unmanageably large. It is necessary to optimize or 'prune' this set down 
to a manageable number [5]. 

The operator set is first pruned based on intrinsic noise. This amounts to selecting the 
ten least noisy operators within each operator type (single-site, singly-displaced, etc.). 
These operators are determined by examining the diagonal correlators (0| Giit) ^(0) |0) . 
We find that these least noisy operators have good overlap with the low-lying states of 
interest and couple weakly to the higher-lying energy modes of the theory. 



In addition, noise enters if the operator set produces a correlation matrix which is ill- 
conditioned. After the least noisy set of operators of each type is determined, the subset 
of those operators with the lowest condition number is chos en. This cond ition number 
is calculated for the normalized correlation matrix Qj{t') / yJCu{t')Cjj{t'), where Qj(t) 
is defined in Eq. 1 and evaluated at some early time t' ~ 3a t . In determining the optimal 
subset of operators, both condition number and operator type diversity are taken into 
consideration. 

An exploratory calculation designed to test the effectiveness of the method has been 
performed for the Nucleon [5] and Delta baryons. To this end, calculations were carried 
out on a small volume (ss 1.2fm) lattice with a relatively coarse O.lfm) lattice spac- 
ing. For computational simplicity, the pion mass for this calculation was unphysically 
large («% 700MeV) and the quenched approximation (omission of quark loops) was 
employed. Also, these results were based on a low- statistics ensemble of 200 configura- 
tions. 

In Fig. 2, the spectrum for each lattice J p value (Gi g ,H g , etc.) is plotted in lattice 
units. Comparison is made with experimental results [6] by identifying these lattice spin- 
parity values with their continuum counterparts using the group theoretical process of 
subduction. It should be noted that for J > |, one continuum state will have several 
lattice counterparts. The scales have been set using the mass of the lowest-lying hadron, 
the Pi i (939). Shown here are the experimental and lattice results for the Deltas. For an 
analogous plot of the Nucleons, see Ref. [7]. 

Upon examination of Fig. 2, we notice that the lowest-lying energy state corresponds 

to a J p = | + (H g ) state, in agreement with the experimental i>33(1232). The next two 

states occur in the j (Gi u ) and | (H u ) channels, in agreement with the experimental 
results . However, the first even parity excitation, the P33 ( 1 600) , disagrees with the exper- 
imental result. This feature is also present in the nucleon sector [7], where the Pi 1 (1440) 
Roper resonance is higher than the first odd-parity excitation. As our simulation ne- 
glected quark loops and possessed an unphysically high pion mass, this discrepancy 
indicates that either these states are not interpolated by standard 3 -quark operators (in 
the quenched approximation) or that they are sensitive to quark mass effects. 

The exploratory calculation described here will be followed by unquenched high- 
statistics runs at larger volumes, finer lattice spacings, and lower quark masses. Com- 
putational time has already been allotted for such runs using 10 million core-hours on 
the Cray XT3 at Oak Ridge National Lab and 3 million service units on the Cray XT3 
'bigben' at the Pittsburgh Computing Center. 

Future work will utilize 'all-to-all' propagators [8], which estimate quark propagators 
from all initial sites to all final sites. This is in contrast with the currently used 'point-to- 
air propagators which only employ one initial site. This method will result in an increase 
in statistics, but more importantly enable the construction of multi-hadron operators. 

In conclusion, this exploratory calculation of the low-lying Nucleon and Delta spectra 
by the LHPC demonstrates the extraction of an unprecedented number of excited states 
in Lattice QCD. With the allocation of future computer resources and the development 
of all-to-all propagators, our collaboration continues to move toward the ultimate goal 
of identifying states with definite quantum numbers, predicting their masses, and calcu- 
lating matrix elements. 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 




Gig Hg G2g 



Glu 



Hu 



G2u 



2.8 - I, ,,(2750) I, n (2750) (4) 

L n (2750) (x2) 
H 3|1 (2420)(x2) 3J3 

H 3.u< 242 °) H 3 |1 (2420) G ^ (2400) 0^(2400) (x2) 

F„(2390) 

_H,.(2300) H,„(2300) (x2) G, 7 (2200) 



2.4- 



D„(2350) 



2 _ F 37 (1950) 
P 3| (1750) 



F„(2000) 



S„(2150) 



> 1.6 
O 

s 

1.2 ■ 



0.8- 



i i D 33 (1940) 

' F~05)s 3 ~ ) D 35( 193 °) 
"^71700) 



P„(1600) 



P„(1232) 



S 31 (1620) 



0.4- 



Glg 



Hg 



G2g 



Glu 



Hu 



G2u 



FIGURE 2. The low-lying Delta spectrum. The left plot is the lattice calculation (colors are used to 
distinguish different levels). In the experimental plot (right), black denotes a 4-star state, blue a 3-star 
state, brown a 2-star state, and gray a 1-star state. States are grouped according to their lattice spin-parity 
values, Gig,H g , etc. 
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